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Abstract 

We study the isochronicity of centers at O G R 2 for systems x = — y + A(x,y), y = x + B(x,y), 

where A, B E R[ac, y], which can be reduced to the Lienard type equation. Using the so-called 

C-afeorithm we have found 27 new multiparameter isochronous centers. 
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1 Introduction 



1.1 Generalities 

Let us consider the system of real differential equations of the form 

= x = -y + A(x,y), -| = y = x + B(x, y), (1) 

where (x,y) belongs to an open connected subset U C R 2 , A,B G C l (U, R), where A and 
B as well as their first derivatives vanish at (0,0). An isolated singular point p G U of 
system (CQ) is a center if there exists a punctured neighborhood V C i7 of p such that every 
orbit of (CQ) lying in V is a closed orbit surrounding p. A center p is isochronous if the 
period is constant for all closed orbits in some neighborhood of p. 

The simplest example is the linear isochronous center at the origin O = (0, 0) given by 
the system 

x = -y, y = x. (2) 

The problem of caracterization of couples (A, B) such that O is an isochronous center 
(even a center) for the system ([T]) is largely open. 

An overview [I] present the basic results concerning the problem of the isochronicity, 
see also [H D2] . 

The hunting of isochronous centers is now a flourishing activity. By this paper we would 
like to contribute to it. 



The well known Poincare Theorem asserts that when A and B are real analytic, a 
center of ([T]) is isochronous if and only if in some real analytic coordinate system it take 
the form of the linear center (j2j). Let us formulate now another theorem of the same vein 
(see for example [lj, Th.13.1 and [17], Th.4.2.1). 

Theorem A (|15j. Th.3.3) Let us suppose that the origin O is an isochronous center of 
system ([T]) with real analytic functions A and B. Let F(x,y) = x 2 + y 2 + o(\(x,y)\ 2 ) be 
an real analytic first integral defined in some neighborhood of O. Then there exists a real 
analytic change of coordinates u(x,y) = x + o(\(x,y)\), v(x,y) — x + o(\(x, y)\) bringing the 
system §T§ to the linear system u = —v, v = u and such that F(x,y) = u 2 (x,y) +v 2 (x,y). 

We now pass to the heart of the matter. To make this paper more accessible, we report 
all strictly technical remarks concerning C-algorithm and Grobner basis to Appendix, 
Sec. 7. 
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In some circumstances system ([I]) can be reduced to the Lienard type equation 

x + f{x)x 2 + g{x) = (3) 

with /, g G C 1 (J, R), where J is some neighborhood of G R and g(0) = 0. If it is 
so, the system (CQ) is called reducible. To the equation (j3J) one associates equivalent, two 
dimentional (planar), Lienard type system 



x = y 
V = -g(x) - f(x)y 2 



(4) 



For reducible systems considered in this paper, the nature of singular point O for both 
system (JT]) and (j3J) is the same; in particular this concerns the centers and isochronous 
centers. More precisely, for the purpose of this paper we shall consider two cases where 
such a reduction is possible. 

• Case 1: When —y + A(x, y) = —yA(x) and x + B(x, y) = B(x) + C(x)y 2 system (TjQ) 
can be written 

x = —yA(x) 
y = B(x) + C(x)y 2 

By the change of coordinates (u,v) := (x, —yA(x)) we get 



(5) 



u = v 



■r , s A'{u) - C(u) , 

v = -A(u)B(u) H —v 2 

A(u) 

In this way we obtain the reduction to the system ([3]) with 

= _ A'(x) -C(x) ^ g{x) = A{x) z {x) _ (6) 
A(x) 

Case 2: When A(x, y) = and B(x, y) = xP(y) where P(0) = 0. In this case system 
(ED can be written 



x = -y 
y = x(l + P{y)) 
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By the change of coordinates (u,v) := (y,x(l + P(y))) we get 



u = v 

P'(u) 



v = -u(l + P{u)) + V 
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1 + P(u) 
We obtain the system (jHJ) with 

f( x "> = ~ i \ and = < l + p ( x ))- ( 7 ) 
1 + P{x) 

In both cases the determinant of the Jacobian matrix of coordinate change does not 
vanish at (0, 0). Thus the nature of singular point O is the same for system (CQ) and (j4]). 
Let us return now to the Lienard type equation ([3]). Let us define the following functions 

PX PX 

F(x) := / f(s)ds, (j)(x) := / e F{s) ds. (8) 
Jo Jo 

When xg{x) > for x ^ 0, define the function X by 

\x{xf= P g(s)e 2F{s) ds (9) 



^ Jo 

and xX(x) > for x ^ 0. 

Let us formulate now the following theorems which are the starting point of this paper. 
Theorem B ([18], Th.l) Let f, g e C\J,R). If xg(x) >0forx^0, then the system @> 
has a center at the origin O. When f, g are real analytic , this condition is also necessary. 

When f,g£ C 1 (J, M), the first integral of the system (j3j) is given by the formula 

I(x, x) = 2 f g(s)e 2F{s) ds + {xe F{x) f (10) 
Jo 



Theorem C([T0], Th.2.1) Let f, g be real analytic functions defined in a neighborhood 
J of G M., and let xg(x) > for x ^ 0. Then system (j4j) has an isochronous center at 
O if and only if there exists an real analytic odd function h which satisfies the following 
conditions 

XU) 9{x)e F ^\ (11) 



1 + h(X(x)) 
the function <f>(x) satisfies 



fX(x) 

0(x) = X(x) + I h(t)dt, (12) 



and X(x)<f)(x) > for x^O. 

In particular, when f and g are odd, O is an isochronous center if and only if g(x) = 
e~ F ^(p(x) , or equivalently h = 0. 

The function h is called Urabe function. The above Theorem implies 
Corollary A ([10j. Corollary 2.4) Let f , g be real analytic functions defined in a neighbor- 
hood o/OgM, and xg(x) > for x ^ 0. The origin O is isochronous center of system (j4]) 
with Urabe function h = if and only if 

g'(x)+g(x)f(x) = l (13) 

for sufficiently small x. 

In the future we shall call the Urabe function of the isochronous center of reducible 
system (CQ) the Urabe function of the corresponding Lienard type equation. 

In [10] the second author described how to use Theorem C to build an algorithm (C- 
algorithm, see Sec. 7.1 Appendix for more details) to look for isochronous centers at the 
origin for reducible system ([T|), and apply to the case where A and B are polynomials of 
degree 3. This work was continued in [TT] . 

The main results obtained in [TO] and [TT] are the necessary and sufficient conditions 
for isochronicity of the center at O in term of parameters for the cubic system 

x = —y + axy + bx 2 y ! 
y = x + a\x 2 + a^y 2 + a 4 x 3 + a^xy 2 I 

The aim of this paper is to extend these investigations for systems with higher order 
perturbations of the linear center x = —y, y = x. 

Like in [Ml EH], our main tool to investigate the isochronous centers for multiparameters 
systems reducible to Lienard type equation is C-algorithm. Nevertheless, when searching 
only the isochronous centers with zero Urabe function the Corollary A gives a much simpler 
method which is widely used in this paper. It consists in identifying the parameters values 
for which identity (TT3"j) is satisfied. 

In all cases considered in [TUl [TT] as well as in the present paper the Urabe function 
is of the form h(X) = /; T ^„ where s is an odd natural number, ki,k2,k 3 G R 
and /c2 > 0. Like in [TUl ED]; we as k if the Urabe function of corresponding Lienard type 
equation (called in the sequel also the Urabe function of the isochronous center under 
consideration) is always of the above form. 
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One of our contributions is the explicit description of simple multiparameter families of 
system (pQ) with isochronous centers at the origin and with a very complicated coefficients. 
Their complexity clearly indicates that we approach the end of purely enumerative study 
in this field. 

Let us stress that using the change of variables given by a polynomial automorphism 
of M. 2 it is easy to transform a simple system of polynomial differential equation with 
isochronous center at the origin into a very complicated one. But systems thus obtained 
do not belong to the class of simple and natural systems studied in the present paper. Our 
contribution is the explicit description of such complicated systems in simple and natural 
multiparameter families of planar polynomial differential systems. 

In our investigations we have used Maple in its version 10. To compute the Grobner 
basis (with DRL order) of the obtained systems of polynomial equations, we have used 
Salsa Software more precisely the implementation FGb \12\ . 

1.2 Beyond the degree 3 

We now present the list of reducible systems for which we study the isochronous centers 
at the origin. 

1. In Section 2 we study the most general homogeneous perturbation of arbitrary degree 
n > 3 of the linear center which belongs to the Case 1 from the Sec. 1.1 : 



Here we found 3 isochronous centers for even n > 4 and 2 isochronous centers for 
odd n > 3 which are new. 

2. In Sections 3 and 4 we study the most general polynomial perturbation of degree 
four of the linear center which belongs to the Case 1 from the Sec. 1.1 : 



First using Corollary A we identify all isochronous centers with zero Urabe function. 
Here we found 6 isochronous centers which are new. The study of this system by 
C-algorithm can not be performed by our actual computer facilities. Thus, we select 




(14) 



x = —y + a\\xy + a 2 ix 2 y + a 3 ix 3 y 

y = x + b 20 x 2 + b 30 x 3 + b 02 y 2 + b 12 xy 2 + b 22 x 2 y 2 + b 40 x 4 



(15) 
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for investigation two sub-families; the first one when ai,i = 63,0 = and the second 
one when a^i = a 2> i = 0. Here we found 10 isochronous centers which are new. 

3. In Section 5 we study the most general polynomial perturbation of degree five of the 
linear center which belongs to the Case 1 from the Sec. 1.1 : 

x = —y + a u xy + a 21 x 2 y + a 3 ix 3 y + a^x y 

y = x + b 20 x 2 + b 30 x 3 + b 02 y 2 + b 12 xy 2 + b 22 x 2 y 2 + b 32 x 3 y 2 + 640a; 4 + b m x^ 

Using Corollary A we identify all isochronous centers with zero Urabe function where 
650 = 0. Here we found 8 isochronous centers which are new. 

4. In Section 6 we study the following Abel system of arbitrary degree n > 2 which 
belongs to the Case 2 (see Sec. 1.1) : 

x = -y 

n 

V = ^2a k xy k , 

k=0 

where a k G M, for k = 0, . . . ,n. Here we verify that up to n = 9 there are no other 
isochronous center than the one found in [19]. 

To sum up, we have found 24 multiparameter isochronous centers as well as three 
infinite families of them that correspond to the perturbations of arbitrary high degree, the 
whole of which are new. 

Concerning the reduction to the Lienard type equations the systems f|T4^ - f|T6|) come 
under case 1, while system ffTTl) come under case 2 (see Sec. 1.1). In particular, for the 
systems f|T4|) - ffTB]) the functions / and g from equation (3) are those given by formulas (jSJ), 
while for the system (|T7|) they are those given by formulas (171) ■ 

Let us stress that by Theorem B, in all the above cases the origin O is always a center 
(indeed, the condition xg(x) > for x 7^ is satisfied for sufficiently small 

When describing in Sec. 3-6 the identified isochronous centers, all parameters intervening 
in the formulas are arbitrary, except that one always supposes that the denominators are 
non zero. To avoid misprints all formulas are written exactly in the form produced by 
Maple. All fractions which appear in the formulas are irreducible. In all cases when we 
were able to write down first integrals and linearizing changes of variables, the explicite 
formulas are reported. 
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2 Homogeneous perturbations of arbitrary degree 

Taking into account the condition g'(x) + f(x)g(x) = 1 from Corollary A, one easily obtains 
the following Theorem 

Theorem 2.1. For n > 2 the system (I14p has an isochronous center at the origin O with 
zero Urabe function only in one of the following two cases 



x = —y + ax n l y 
y = x + ax n ~ 2 y 2 



x = y + -x n l y 



n 

y = x + bx n ~ 2 y 2 — ^ n ~P b x n 
Moreover, for odd n > 3 t/iere are no other isochronous centers. 

Proof. System ffT4l) is reducible to system (j4]) with 

= x"- 2 (6 + na -a) ^ = ^ _ ^ 
v ' 1 - ax™" 1 v 7 v ' y ' 

The condition g'(x) + f(x)g(x) = 1 allows directly to the following two cases : 

1. {a = b, c = 0} which gives the system (fT8|) . 

2. -fc = — O = ^ 1 which gives the system ( fl9l) . 

Applying formula ffTU]) using Maple, one see that for n = 4, 6, 8 the first integral of sys- 
tem ffl8|) takes the form 

_ (x 2 + y 2 ) 

n Si8\ — [XT 

(-1 + ax™" 1 )™- 1 

Then Theorem A suggest that the linearizing change of coordinates is 

u = u X v = - V (20) 



(18) 



(19) 



\/l - ax"" 1 ' u -\/r 



a x 



n-l 



Now one directly verifies that is always a first integral of system ( Ti~8i) and using Maple 
one easily checks that (120!) is a linearizing change of coordinates. 

Exactly the same arguments work for the system (fl9|) . Its first integral is 

x 2 (1 + ex"" 1 ) 2 + y 2 



2 2 71 

n — 1) (n — 1 + na"^ 1 )"- 1 



and its linearizing change of coordinates is 

x (1 + cx n ~ l ) 



u 



(n — 1) ((n — 1 + ncx n (n — 1) ((n — 1 + ncx n 



□ 



When n > 4 is even the preliminary investigation of system (1141) performed by C- 
algorithm strongly suggests that for such n there exists exactly one additional isochronous 
center with non zero Urabe function. Its existence is proved in Theorem l2.21 Unfortunately, 
its uniqueness is not yet proved for arbitrary even n > 4. For n = 4, 6, 8 the uniqueness 
was proved using Maple and Grobner Basis method. 

Let us point out that our final proofs are done by hand computations, without using 
computer algebra. 

Theorem 2.2. The system ffT41) with arbitrary even n > 2 and a = 2b, c = —b, b ^ 0, 

has an isochronous center at the origin with non zero Urabe function 

h(X) bX "- 1 



VI + b 2 X 2n ~ 2 

Proof. When a = 2b and c = —b the system (Tl4l) becomes 



x = —y + 2bx n x y 
y = x + bx n ~ 2 y 2 — bx n 

The change of variables (x, y) \ — > (x/b, y/b) reduces the system ( l2lT) to the form 

x = —y + 2x n ~ 1 y 
y = x + x n ~ 2 y 2 — x r ' 

which is reducible to the Lienard type equation ([3]) with 



(21) 



( — 1 4- 2 n ) r n ~ 2 
f(x) = 1 1 _ 2x L and g(x) = (l - 2 x^ 1 ) (x - x n ) 



Then 



F(x)= / ( s )d s = __-ln(l-2x"- 1 ) 
which gives the right hand side of the equality (TTTT) 

g(x)e^ = X{1 X ^ \ . 



On the other hand, e 2F( ^ = (1 — 2x n x ) . 
^From the equation (jUJ) we compute 



x 



X(x) = J 2 J q 9(s)e 2F ^ds 



n-l^ 



o (1 — 2 a; 
and 

X(x) n ~ l x n ~ l 

h(X(x)) = 



1 + X(x) 

Then we compute the left hand side of the equality ([TT|) : 



X(x) (l- 2 j)"- 1 )^2 x(1 - x n x ) 



l + h{X(x)) l + j^r ~ (l-2x«- 1 )^ 

Which proves that the equality (TTTT) is satisfied. Let us stress that the above computations 
remain valid for every n > 2. Nevertheless, for n odd /i is not an odd function and thus it 
is not an Urabe function which is odd by definition. □ 

Theorem 2.3. For arbitrary n>2, the system (I2TI) has the following first integral 

Hm - 2bx^ - 1 • 

Proof. Using formula (flOl) . one easily computes by Maple the first integral for n = 4, 6, 8. 
The obtained results strongly suggest the veracity of the formula for -f/jgrj. Now one easily 
can check by hand that H^n is a first integral. □ 

Let us return to system ()14p . It is well known that for n = 2, this system has an 
isochronous center in exactly four cases, so called Loud isochronous centers (see [Til [TU] ) . 
They correspond to (a = b, c = 0), (a = ~, c = — |), (a = 26, c = —6) and (6 = |, c = 0). 
The first two are those from Theorem 12.11 , the third is the one from Theorem 12.21 

Let us note the Taylor expansion of the Urabe function h(X) = C\X + C3X 3 + . . .. 
As noted at the begining of the Section, for n = 4 one has exactly 3 cases of isochronous 
centers. Why such a difference? The difference is in the algebraic structure of the equations 
generated by C-algorithm. For n = 2, the second of such equations is — 3ci + a — 2c — b = 
and Ci can be non zero, while for n > 3, the second such equation is always c\ = 0. Thus 
the freedom for existence of non zero Urabe function is greater for n = 2 than for n > 3. 
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3 Non-homogeneous perturbations of degree four with 
zero Urabe function 



Taking into account the condition g'{x) + f(x)g(x) = 1 from Corollary A, using Maple one 
easily obtains the following Theorem. 

Theorem 3.1. The system ( fl5l) has an isochronous center at the origin O with zero Urabe 
function only in one of the following six cases, where one supposes that all denominators 
are non zero polynomials. 



x = -y + b 02 xy + a 2l x 2 y + a 31 x 3 y 
y = x + b 02 y 2 + a 2 ixy 2 + a 31 x 2 y 2 



II 



x =-y + b Q2 xy - § b 30 x 2 y + b 02 b 30 x 3 y 

y = x + b 02 y 2 + b 30 x 3 - § b 30 xy 2 + 3 b 02 b 30 x 2 y 2 



III 



x = -y + (&02 + 2 6 2 o) xy + a 21 x 2 y + (b 20 a 21 - b 02 b 20 2 - 4 6 20 3 ) x 3 y 

y = x + b 20 x 2 + b 02 y 2 + (a 21 + b 02 b 20 + 4 b 20 2 ) xy 2 
+ (2 b 20 a 21 - 2 b 02 b 20 2 - 8 b 20 3 ) x 2 y 2 



IV 



(-9f)30&20 2 -b2() 2 a2i+26 2 o 4 +4f)3oa2i+6fe30 2 ) 

x = —y + - -, 5 1 L xy 

y 620 (-620 2 +4 630) ^ 



9 63o(-2b3ofe20 2 -t | 20 2 a2i+4b3oa2i+663o 2 ) r, 

+a2ix y + — — ^(-^o 2 ^^) ~ x y 

7 2 (-17b30fr20 2 -620 2 a2i+4&20 4 +463 a 2 i+6f>30 2 ) 2 

OonX H - — 7 — ■ — o — - — x y 



y = x + b 20 x z + 



&2o(-&20 2 +4fc3o) 



j, 3 1 O (- 6 20 2 a2i+4t>3oa2l+b3ob20 2 -3 630 2 ) 2 
+ t)30X +2 -620 2 +4 6 3 ^ 



b3o(-2fe3,ofe20 2 -teo 2 a2i+463oa 2 i+663o 2 ) 9 9 

+3 — * t 5 > L x z y z 

b2o(-b20 2 +4 6 3 o) y 
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V 



(-32 6 4 of>3ob20 2 +42 640^20+8 b 40 b 2 Q 4 +b 4 ob 3 o 2 -2 b 3 2 & 2 3 +7 &30 3 fc 2 o) 
X ~~ _ ^ -fe 3 o 2 feo 2 +4 a 2 , 4 ,o620 3 -18 640630^0+27 6 4 o 2 +4 630 3 



(6 6 3 o 4 -2 f>20 2 &30 3 -27 f)4oteofe30 2 +8 fe 3 ob4ofe20 3 +39 f)3ofe 4 o 2 -4 fe 40 2 fe20 2 ) 2 
- 630 2 &20 2 +4 640 &20 3 - 18 640 630 &20 +27 640 2 +4 fe 3 3 X y 



„ &40(-b 3 2 620 2 -14fe40''30f'20 + 18t'40 2 +4b40fe20 3 +3fe30 3 ) 3 
~~ - ho 2 620 2 +4 640 feo 3 - 18 640 &30 &20 +27 640 2 +4 b 30 3 ^ ^ 



, 2 (-68 6 4 0b30fe0 2 +96b40 2 620 + 16 640fc20 4 +'>40fe30 2 -4fe30 2 fe20 3 + 15 630 3 b20) 2 

X + 2 ° X - 630 2 &20 2 +4 fe 40 620 3 - 18 b 40 630 620 +27 6 40 2 +4 630 3 ^ 



+6 30 ^ 3 - 2 



(9 b 3 Q 4 - 3 b 2 Q 2 b30 3 - 40 fe 40 &20 b 3 Q 2 + 1 2 b 3 Q ^40 fe 2 Q 3 +60 b 3 Q &4Q 2 - 8 fe 4 Q 2 fe 2 Q 2 ) 2 
-b 3 O 2 b2O 2 +4fe4ob20 3 -18fe4o630b20+27b40 2 +4 630 3 ^ 



b40(-b30 2 fe0 2 -14fe4of>30f)20 + 18f)40 2 +4640fe20 3 +3fe30 3 ) 2 2 , t 4 
-f*3O 2 &20 2 +4b40fc2O 3 -18640630b20+27fc 4 O 2 +4b3O 3 ^ ^ + ^ 



VI 



X 



y 



= -y + (2 b 20 + 602) xy - ^ (567 Z 2 6 20 + 24 Z6 02 - 804 Z6 20 + H3 6 20 - 8 b 02 ) x 3 y 

feo (13032 Z 2 fe 2-4488 Zb 2+68719 Z 2 6 20 -22970 Zb 20 +384 6 2 +1943 b 20 ) 2 
^ 24(387 Z 2 -144 Z +13) X V 

._,,), T 2 , l 2 Z6 20 2 (27Z-17) 3 
— X -+■ O20X -+■ Oo2y 4+12Z X 

- b -§ (-108 Z 2 b 02 + 1053 Z 2 6 20 + 152 Zb Q2 - 1964 Z6 20 - 76 6 02 + 155 6 20 ) xy 2 



-bm_ ( 567 Z 2 &2o + 24 ^ _ g04 zh ^ + 113 ^ _ g 6o2 ) ^2 + 1/4 ^^4 
where Z is the only real root of the equation 27 s 3 — 47 s 2 + 13 s — 1 = 0, which is equal to 
Z 



1 



+ 



47 



1 3/ 7= 1156 
— V 39428 + 324 \/93 + = 

81 81 V39428 + 324 V93 81 



4 Non-homogeneous perturbations of degree four 

Let us consider system (fT5j) . We would like to identify all its isochronous centers by C- 
algorithm, without taking into account the nature of its Urabe function. In full generality, 
this problem cannot be attained by our actual computer possiblities. Indeed, we do not 
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succeed to compute a Grobner basis for the nine C-algorithm generated polynomials on 9 
unknown {ay} and {b sr } of all even degrees between 2 and 18. 

Inspecting the system under consideration one sees that the annulation of some param- 
eters {ciij} and {b rs } will substantially simplify the system. This is the reason of our choice 
of two families presented below. 

4.1 First family 

Let us assume an = b 30 = 0, in this case system (fT5l) reduces to the system 



x 

y 



-y + a2ix 2 y + a 31 x 3 y 



x + b 20 x + b 02 y + bnxy + b 22 x y + b i0 x 



(22) 



Theorem 4.1. The system (|22|) has an isochronous center atO if and only if its parameters 
satisfy one of the folowing 6 conditions : 





x 

y 



y + b -f x 3 y 



x + b 22 x 2 y 2 - x 4 



x = —y + 2b 22 x 3 y 

y = x + b 22 x 2 y 2 - b 22 x 4 



x = —y + a 2 \x 2 y + a 3 ix 3 y 
y = x + a 2 \xy 2 + a 3 ix 2 y 2 



IV 



x 

y 



-V+'i b 20 2 x 2 y - I b 2 o i x 3 y 
x + b 20 x 2 - 2 b 20 y 2 + § b 20 2 xy 2 - | b 20 6 x 2 y 



3^2„,2 



V 



x = -y + § b 20 2 x 2 y - § b 20 3 x 3 y 
y = x + b 20 x 2 - 2 b 20 y 2 + § b 20 2 xy 2 



§ b 20 3 x 2 y 2 



VI 
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x = -y + b 20 2 (2 + gfe) x 2 y + a 31 x 3 y 1 
y = x + b 20 x 2 - 2 b 20 y 2 + b 20 2 (4 + g^) xy 2 + 2 a 31 x 2 y 2 J 
where b 2 Q 7^ 0. 

Proof. C-algorithm gives the six candidates to be isochronous centers. We had to derive 
19 times to get the necessary conditions of isochronicity. 

To apply succesfully the C-algorithm, we use the two tricks explained in Appendix, 
Sec. 7.2 : homogenization and reduction of the dimension of the parameters space by one. 
This leads to the proof that the cases I- VI of Theorem 14.11 satisfy the necessary conditions 
of isochronicity. We check that the necessary conditions are also sufficient by direct applica- 
tion of Corollary A to the cases I, III- VI. Indeed, in all those four cases g'(x)+f(x)g(x) = 1. 
For sufficiently small x the case II is a particular case of the system (12"TI) when n = 4, studied 
in Theorem 12.21 □ 

Let us note that among the above six cases only the cases I, II and III with 021 = 
represent the homogeneous perturbations. All other cases are non-homogeneous. 

Note also that the above three homogeneous cases were already identified by Theorems 
12.21 and 12.11 But contrary to the quoted Theorems, here we have the exhaustive list of 
isochronous centers for n = 4. 



4.2 Second family 

Consider system f)15p . with an = a 2 i = 0. We obtain the seven parameter real system of 
degree 4. 



x = -y + a 31 x 3 y 



y = x + b 20 x z + b 02 y z + b 30 x 3 + b l2 xy z + b 22 x 2 y 2 + b 40 x 



(23) 



Theorem 4.2. The system (123]) has an isochronous center atO if and only if its parameters 
satisfy one of the folowing seven cases : 

The three cases I, II and III with a 2 i = come from Theorem \4-l\ and correspond to 
homogeneous perturbations. The following four cases correspond to the non-homogeneous 
perturbations. 
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X 

y 



-y + j bo2 3 x 3 y 



X 



I b 2X 2 + b 02 y 2 + \ b 02 2 xy 2 + \ b 02 3 x 2 y 2 



V 



x 

y 



= ~y + lk b M (-21 + 5V33)ar 8 y 

= x-\ b 02 x 2 + b 02 y 2 + i b 02 2 (9 - V33) x 3 

+± b 02 2 (-1 + v/33) xy 2 + i 6 02 3 (-21 + 5 v^3) xV 



VI 



± =-y~^ 2 b 02 3 (21 + 5V33)x 3 y 
y = x - | 6 02 x 2 + W + &02 2 (9 + VS3) 

-i 6 02 2 (1 + \/33) xy 2 - i b 02 3 (21 + 5 v/33) xV 



VII 



2 fc 20 3 (-43t 2 / 3 -7670 V3297+12112 ^i+52 ^IV3297-336886) 3 

x = — y + 3549 ^273 ^ 1/ 

_ , 9 1 62 of-3822 62 ot 2/3 -6242964fe 2 0-127764 b 2 ov / 3297+159432b 2 ^t) o 

y = x - 2 6 2> or - Yog^ — * ■ Wa L xy 

1 62, o(l032 6 2 ,o 2 t 2/3 +184080 6 2 ,o 2 v / 3297-290688 f> 2 , 2 v^-1248 b 2 ,oV3297 $^+8085264 6 2 ,o 2 ) 2 2 \ 

— — ^2/3 ^ 1/ ' 



10647 



, 2 j b 2 ,o (-53144 6 2 ,o v^+2080988 6 2 , +42588 62,0 V3297-5824 b 2 ,o< 2/3 ) 3 

+ b 2,oX - Y0647 7573 



-X 



l 6 2 ,o(-2150b 2 , 2 t 2/3 -11926fe2,o 2 ^7+234 6 2 ,o 2 V3297 v / t+1085248 b 2 ,o 2 +18720 b 2 ,o 2 V3297) 4 
10647 W/3 X 



where t = 22868 + 468 V3297 

Proof. Thanks to C-algorithm we obtain the necessary conditions for the isochronicity of 
the center at the origin for system (1231) and we establish the seven cases given in the 
theorem. 

We check that the obtained necessary conditions are also sufficient by direct application 
of Corollary A to the cases IV- VII. Indeed, in all those four cases g\x) + f(x)g(x) = 1 for 
sufficiently small x. 

□ 



The centers I — III of Theorem 2 have been already identified in [TJ. 
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5 Non-homogeneous perturbations of degree five with 
zero Urabe function 

By Corollary A the problem is reduced to solving the equation g'(x) + f(x)g(x) = 1, with 
/ and g defined in Sec. 1.1 with respect to the system (1161) . In this case the equation 
g'(x) + f(x)g(x) = 1 is equivalent to some system of 8 polynomials depending on 12 
unknown {a^} and {b sr } of degree 2, 2, 2, 1, 2, 2, 2, 2. Applying the Grobner basis method 
one obtains a basis of 90 polynomials whose degrees varies between 1 and 8. This system 
is too hard to handle. By inspecting the system one sees that when 6 50 = it is very much 
simplified. This is the reason of our choice 650 = 0. 

Theorem 5.1. The system ( fl6l) where 650 = has an isochronous center at the origin O 
with zero Urabe function only in one of the following 8 cases, where one supposes that all 
denominators are non zero polynomials. 



x = -y + a n xy + b 12 x 2 y + a 3 ix 3 y + b 32 x 4 y 1 
y = x + any 2 + bnxy 2 + a 31 x 2 y 2 + b 32 x 3 y 2 J 

77 

x = -y + a n xy + a 31 x 3 y - | a n a 31 x 4 y 1 
y = x + a 1± y 2 + 4 a 31 x 2 y 2 - § a 3Y x A - 3 aua 31 x 3 y 2 J 

777 

x =-y + f^xy + (3b 30 + b 12 )x 2 y + a 31 x 3 y+(%b 30 2 + b 12 b 30 )x*y 1 
V = x + b^y 2 + & 3ox 3 + b l2 xy 2 + 3 a 31 x 2 y 2 + (f b 30 2 + 3 612630) x 3 y 2 J 



IV 

x = -y + a u xy + (612 - 2 6 2 o 2 - 011620) z 2 2/ + a 3 ix 3 y+ 
(-612620 2 + 4 620 4 + 2 an6 2 o 3 + 031620) x 4 y 

> 

y =x + b 20 x 2 + (-2 6 20 + a u ) y 2 + b 12 xy 2 + (6i 2 6 20 - 4 6 20 3 - 2 a u b 20 2 + a 31 ) x 2 y 2 
+ (-2 6 12 6 20 2 + 8 6 20 4 + 4 a n 6 20 3 + 2 a 3 i6 20 ) x 3 y 2 
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V 



x = -y + a n xy + a 31 x 3 y 



VI 



(l3(>3o 2 b20-ll 6 3o620 3 -5b3oailb20 2 +2fe20 5 +aii620 4 +«3l^20 2 -4a3lb30+4aiife30 2 ) 2 

7 2^ ^ y 

b'20 ^4 bso — 620 J 

&30( — t | 30 a ll&20 2 +7&30 2 &20 — 2 ^O^O 3 — 4 131 ^30+4 ail b30 2 + a 31 ^20 2 ) A 

f>2o (4^)30— teo 2 ) 



y = x + &2o^ 2 + (-2 620 + an) + & 3 o^ 3 



(25 6 3 o 2 620-22f)3ofe20 3 -9fe3oailb20 2 +4b20 5 +2aiib20 4 +a3ib20 2 -4a3ife3o+4aii63o 2 ) 2 

62o(4 6 3 o-02o 2 ) ~~ ''^ 

(763o 2 020-263ob20 3 -b3oaiifeo 2 -2a3i620 2 +8a3i03o+4an63o 2 ) 2 2 
4fe30-»20 2 ^ 

g fe3o(— &3oailb20 2 +7b30 2 f'20 — 2b3oteo 3 — 4a3ife3o+4an63o 2 +a3ife20 2 ) 3 2 
f>2()(4o30-f>20 2 ) ^ ^ 



(108 fc 40 2 -42 f> 40 &20 3 +81 a 3 lb40+teo 6 -3o20 3 a3l) 
^ ~ ~ V + fe20 2 (-&20 3 +27640) 



(-36 4 o&20 3 -fe20 3 a3i+27a3ife 4 o+36fe40 2 ) 9 q 

+ 3 b 2o(- fe 2o 3 + 2 7b40 ) ~ x y + a3i;r v 



q fe 4o(-020 3 a3l+36fe 4 o 2 +27a3ife 4 o) 4 
+ 620 2 (-6 2 o 3 +27b 40 ) ^ 27 



, 2 q ( &20 6 - 32 b 40 bzo 3 +36 b 4 u 2 +27 a 3 i b 40 - b 20 3 a 3 i ) 2 
y -X + O^ +6 b20 2(_ fe20 3 + 27 640 ) ^ 



1 , 2 3 I r> (-4fe 40 b20 3 -b20 3 a3l+27a 3 ib 40 +36 b 4 o 2 ) 2 
+ 3^0* +6 b2() (_ b20 3 +27540 ) X V 



„ (b20 3 a3i-27a 3 ib 4 o+12b 4 o 2 ) 2 2 
_d -feo 3 +27fe 40 X ^ 



7 4 1 1 o 6 4o(-620 3 a3l+36 6 4 o 2 +27a3ife 40 ) g 2 
+640;r + 12 te o 2 (-feo 3 +27b 40 ) ^ ^ 
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VII 



X =-y + anX y-3 2 V 7bV+4&30 3 L Xy 



( 5 630 3 640 +36 640 3 - 2 7 a 1 1 6 30 6 4 o 2 - 4 a 1 1 6 30 4 ) a 
27 640 2 +4 630 3 * 



640 ( 640 630 2 + 2 7 a 1 1 640 2 +4 a 1 1 6 30 3 ) 4 



27 640 2 +4 630 3 



xry 



O , , Q _ 630(3630+20640 ) O 

1/ = x + a n y z + b 30 x 6 - 6 ■ 



27 640 2 +4 630 



(-27aii&3o&4o 2 -4a 11 63(, 4 +7630 3 640+486 4 o 3 ) 2 2 
' d 27 6 40 2 +4 6 3 3 X V 



,7 4 I /I ^4o(64o630 2 +27aiife4o 2 +4an63o 3 ) g 2 
+&40X + 4 27&40 2 +4 6 3 o 3 X V 



VIII 



x = -y + anxy + a 2 ix 2 y + a 3 ix 3 y + a 41 x 4 y 



y = x + b 20 x z + b 30 x 6 + b 02 y z + b 12 xy z + b 22 x 2 y 2 + h 2 x A y z + 6 40 a; 4 + b 50 x 



„32 



where 



Pli _ b 2 oP 2 \ _ P31 _ P41 _ P02 , _ P22 , 

au ~ g"' a21 " ~Q~' a31 " q" ' a41 " q" ' 602 " ~Q ' 622 " g"' 632 
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P u = -41 b 40 b 30 2 b 20 + 2 b 40 b 12 b 20 3 - 9 b 40 b 30 b 12 b 20 + 9 6 30 4 - 8 b 40 b 20 5 
-50 6 4 o 2 & 20 2 + f b 12 b m 2 + 60 6 40 2 &30 - 10 h % 2 + 2 & 30 2 & 20 4 
-\ b 3 o 2 b 2 o 2 bi2 + 44 b 40 b 30 b 20 3 + 2 &i 2 6 3 o 3 

P 21 = -13 b 40 b3o 2 b 2 o + 2 640&12&20 3 - 9 b 4Q b 3Q b 12 b 20 + 3 6 30 4 - 4 6 4 o 2 ^2o 2 

+ f &12&40 2 + 21 b 40 %o - b 30 3 b 20 2 - \ b 30 2 b 20 2 b 12 + 4 b 40 b 30 b 20 3 + 2 &i 2 6 30 3 

P31 = 9 b 30 5 + 2 b 30 % 2 - 3 6 3 oX 2 - 46 b 40 b 20 b 30 3 - | b 30 \ 2 b 20 2 

+ 14 6 3 o 2 &20 3 ^40 ~~ 9 6 3 o 2 &40^12^20 + 60 6 30 2 6 4 o 2 + 2 6 3 o& 40^12^20 

+f &30&12&40 2 + 20 b 30 b 20 2 b 40 2 - 36 6 40 3 &2o - 8 b 20 % 2 
Pa = f &12&40 3 + 60 & 40 3 &30 + 2 b 40 b l2 b 30 3 + 12 b 40 2 b 30 b 20 3 - 8 6 4 oX 2 

I 2 2 2 423 32 

— 9 b 4 ob 3 Q b 2 o b 42 — 9640 b 3 obi 2 b 2 o + 9 640&30 + 2 6 40 z 'bi 2 b 2 Q — 3 6 4 o&30 ^20 

-40 6 40 2 &30 2 &20 

-P02 = -16 b 40 b 20 5 + 80 640&30&20 3 - 104 b 4 o 2 b 20 2 + 4 6 30 2 &2o 4 - 18 &30 3 ^2o 2 
—41 b 40 b 30 b 20 + 2 6 4 o&i2^2o' — 9 b 40 b 30 bi 2 b 20 + 9 6 30 4 + f &12&40 
+60 6 40 2 &30 - \ b3o 2 b 20 2 b 12 + 2 b 12 b 30 3 

P 22 = 180 6 30 2 6 4 o 2 - 144 b 40 % - 32 b 20 % 2 + 6 & 30 &40&i2&20 3 + 27 6 30 5 + 6 6 30 % 2 
-9 b 30 % 20 2 + f b 30 b 12 b 40 2 - 144 b 40 b 20 b 30 3 + 44 b 30 2 b 20 3 b 40 

2 2 2 3 3 2 

+88 6 3 o&2o ^40 — 276 3 o ^40^12^20 — 2 b 3 o b 42 b 2 o 

P32 = 54 b 12 b 40 3 + 240 6 40 3 6 30 + 8 b 40 b 12 b 30 3 + 48 b 40 2 b 30 b 20 3 - 32 6 40 3 6 20 2 
-2 &4ofr30 2 fr2o 2 6i2 - 36 b 40 2 b 30 bi 2 b 20 + 36 b m b m A + 8 b 40 2 b 12 b 20 3 
-I2b 40 b 30 % 2 - lQ0b 40 2 b 30 2 b 20 

Q = 4 640&20 4 - 18 b 40 b 30 b 20 2 + 27 640^20 - & 3 o 2 &2o 3 + 4 6 2 o&30 3 
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6 Abel polynomial system 



By planar Abel system of order n we mean the system 



x = -y 



y = J2 Pk ^ yk 



(24) 



k=0 ) 

where {Pk(x)}o<k<n are smooth functions. 

This section is concerned by the following Abel system 



x=-y 

y = x(l + P(y)), 



(25) 



with P{y) = a\y + a 2 y 2 + a 3 y 3 + .... + a n y n , G M, for k — 0, . . . , n. This is a particular 
Abel system ( 12~4"1) where Pk{x) := a^x, < k < n. 

6.1 Characterization of isochronous centers 

System (1251) is reducible (see Sec. 1.1, Case 2) to the Lienard type equation ([31) with f(x) 
and g(x) defined by ([7]). Definitions ([8]), (Q and Theorems B and C from Sec. 1.1 remain 
valid. Applied to the Abel system (125]) they give : 

Theorem 6.1. The origin O is a center for the system (T25|) . 

The center at O, is isochronous if and only if there exists an odd function h defined in 
some neighborhood o/0 6l which satisfies the following conditions 



and X(x)<f)(x) > for x^O. 

In particular, when P is an even polynomial then the origin O is an isochronous center 
if and only if P = 0. 



X 



1 + h(X) 
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Proof. xg(x) = x 2 (l + P(x)) > for x^O and \x\ small enougth. Then Theorem B implies 
that the origin O is a center of the system (125]) . 
Now 

F(x)= I f{s)ds = -ln(l + P(»), 



X 

■h) 

thus 

4>{ x ) = f e F ^ds = f - dS n ,^ (26) 
Jo Jo 



1 + P{s) 



Then we obtain 

g(x)e F{x) = x{\ + P(x))e- Hl+P{x)) = x 
Following Theorem C, one obtains 

X(x) 
1 + h(X{x)) ~ X 

as well as the identity ( IT2l . 

For the particular case where P is even, it is easy to see that / and g are odd. Theorem C 
thus implies h = and consequently X(x) = x. From ffl2l one deduces that <f>(x) = X(x). 
Then $8% implies that P = 0. □ 

The following paragraph is devoted to illustrate the last theorem by example. 
6.2 An application 

Let us consider the Abel system fl25l) with n = 9 : 



(27) 




with a?; G R, 1 < k < 9. As follows from Theorem I6.1[ the origin O is always a center 
for (ETD. 



Theorem 6.2. T/ie system (1271) /ias an isochronous center at the origin only in the case 

(28) 



x = -y 



y = x + axy + —xy + — xy 



where a e 
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(29) 



Proof. We apply C-algorithm for 

= - 1 s and ^ x ) = < i + p ^)' 

1 + -r (x) 

where P(x) = Yl^=i a i x% - We obtain the unique one-parameter family fl28l) . and computa- 
tions give the Urabe function 

h(X) = = hX 

3 vV + fc 3 X2 

with fci = — ai/3, /?2 = 1, ^3 = 0. 

□ 

By the evident rescalling |x t— > x and fy^y system (BHD takes the form which is a 
particular case of system fl25l) : 

x = -y 
y = x(l + ?/) 3 . 

The isochronous center at the origin O for system (1291) was already depicted in [19] by 
showing that system f[2"9"|) commutes with some transversal polynomial system, but neither 
its first integral nor the linearizing change of coordinates were provided. We shall now 
compute both of them. 

• First integral In the variables u = y and v = x(l + y) 3 the system (1291) is reducible 
to the Lienard type equation u + f(u)u + g{u) = where 

3 

f(u) = and g(u) = u(l + u) 3 . 

1 + u 

By formula (TlTJT) from Theorem B one easily obtains that I(u,v) = n+^p + (i+ u ) a * s 
a first integral of the corresponding planar system 

u = v 

v = -g(u) - f(u)v 2 

Returning to the variables (x, y) one recovers the first integral of the system (J29 

imfav) = %2+ (l + y)2 - 
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Linearization For this purpose we use the method of [13] based on the exitence of 
vector field Y transversal to the vector field defined by the system (1291) and commut- 
ing with it. Maple computations give such a field Y: 



x = x + xy 
y = —x 2 y 3 + y 2 — 3 x 2 y 2 + y — 3 x 2 y — x 2 



(30) 



Following the method described in [13] , we first establish an inverse integrating factor 
V(x, y) of the system ( |29l) : 



V(x, y) = -(y+l) [x 2 + 2x 2 y + x 2 y 2 + y 2 ) 

which leads to the first integrals H (already known) and I of the systems ( 129|) and 
( |30l) respectively: 



H(x, y) = x 2 + 



[i + y y 



I(x, y) = —x + arctan 
Let us define f(z) = z and g(z) = tan(z). 



By the Theorem 4 of [13] we obtain the linearizing change of coordinates 



f(H(x,y))g(I(x,y)) 
v /l + ~g 2 (I(x,y)) 



f(H(x,y)) 



Maple produces the following more explicit formulas that we, as usual, reproduce 
without any change to avoid the misprints: 

" ' \/ T- < « Al \ ■' ^ \ 

u(x,y) 



I -h { I an { x — arctan y ( y+ ^> x 



v(x,y) 



I x 2 + 2 yx 2 + y 2 x 2 + y 



1 + I tan ( x — arctan 



(y+l)x 
y 
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The fact that this change of variables actually is a linearizing one can easily be verified 
by Maple which gives u = — v, i) = u a.s expected. 



Under the light of Theorem 16.21 it is natural to ask if the system ( I28jl is the unique 
system with isochronous center at the origin in family ( 1251) . Even for n — 10, our actual 
computer possibilies are not sufficient to give an answer. 



7 Appendix 
7.1 C- Algorithm 

Theorem C (see Sec. 1.1) leads to an algorithm, first introduced in [TU] (see also [1 lj ) . here- 
after called C-algorithm, which gives necessary conditions for isochronicity of the center at 
the origin O for equation 

Below we recall basic steps of the algorithm. 

Let h be the Urabe function defined in the Theorem C, and u = <p(x). The function 
is invert ible around . 

~ 9{u) ■= TTK(xy (31) 

where now X is considered as a function of u. Our further assumption is that functions 
f(x) and g(x) depend polynomially on certain parameters a := (cti, . . . , a p ) G MP. 

By Theorem C, if the system (|3]) has isochronous center at the origin O, then the Urabe 
function h must be odd, so we have 

oo 

h{X) = Y,C2k + iX 2k+1 

and moreover, 

g{u) = g(x)e F{x) , where x = (p~\u). (32) 

Hence, the right hand sides of fl3Tl) and fl32l) must be equal. We expand both right hand 
sides into the Taylor series around and equate the corresponding coefficients. To this end 
we need to calculate k-th derivatives of f[3T|) and (13"2"|) . 
For (|3~T|) . by straightforward differentiation, we have 

d k g(u) d /d fe -^(w)\ dX 



du k dX V du k ~ l I du 
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Using induction, one can show that for (1321) one obtains 

du k 

where Sk{x) is a function of f(x),g(x) and their derivatives. 

Therefore to compute the first m conditions for isochronicity of system we proceed 
as follows. 

1. We fix m and write 

m 

h{X) = c 2 k-iX 2k ~ 1 + 0(X 2m ), c := (Ci, C 3 , • • • , C 2m _!). 
fe=l 

2. Next, we compute 

do 

Vk ■= W k = S k{0) 

for k = 1, . . . , 2m + 1. Note that those quantities are polynomials in a and c. 

3. By Theorem C we obtain the equations Vk = Wk for k — 1, . . . , 2m + 1. Let us note 
that always Vi = W\ = 1 and thus the first equation is meaningless. 

It appears that we always can eliminate parameters c from these equations. For 
every k > 0, C2k+i occours for the first time, and in a linear way, in the equation 
V2k+2 = W2k+2- This leads to the formula C2k+i = ¥>2k+i(®) for some multivariate 
polynomial </?2fc+i- This brings us in a natural way to the consecutive elimination 
of Ci, c 3 , . . . , c 2m -i- Finally, we obtain at most m polynomial equations si = s 2 = 
S3 = . . . = sm = with p unknowns aij. These equations denoted Sys(m) give M 
necessary conditions for isochronicity of system Q; Sys(m) C Sys(m +1). 

For more details see [TQl [TT] . 

The progammation of the C-algorithm can be done in different ways. Some of them 
appear more successful than others. For the purposes of the present paper we used the 
programm from [2]. 
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7.2 Homogeneization and reduction 

The Grobner bases method for solving the systems of polynomial equations is particularly 
efficient when all these polynomials are homogeneous or weighted-homogeneous (see [3], 



A long experience with C-algorithm indicates that the following facts are always verified 
although they are not proved. 

Let us consider the system (jSJ) written explicitly as 



which as explained in Introduction, is reducible to the equation ([3]) 

When considering the particular cases of system (|33|) corresponding to homogeneous 
(resp. non-homogeneous) perturbations of system (121), the polynomials from Sys(m) are 
homogeneous (resp. non-homogeneous). 

We note that for n — 3, C-algorithm succeeds in establishing isochronicity criteria, 
however for n = 4 the obtained polynomials from the algorithm are much more involved. 
For instance, system fl3"3"|) with n = 4 reduces to the system ffToT) . Now, the first two 
non-zero polynomials from Sys(9) are 



P 2 = 3 a 21 - 3 &i2 + a u 2 - b 20 an - 9 o 30 + 4 b 02 - 5 a u 6o2 + 10 b 2 o 2 + 10 b 20 b 02 , (34) 

P 3 = 72 a 21 2 + 396 b 20 anbi 2 + 90 an& 2&i2 + 36 a lx b 22 + 324 a 3 i6 2 

- 36 a 21 b 12 - 468 6200-11021 + 612 &20O21&02 - 4116 a u b 20 2 b 02 

+ IO8&20O31 - 540 6 30 a 2 i - 3246 40 aii + 1566 6 30 On6 2 - 288 620^22 

- 459 6 30 an 2 - 1296 b 40 b 02 - 306 a 21 a u b 02 + U28b 20 a u 2 b 02 

+ 153a 2 ia,ii 2 - 117 a lx 2 b 12 - 1916 20 oii 3 + 180 ^20^02^12 + 43 an 

- 2319 b 20 a u b 02 2 - 289 a n \ 2 - 360 ^22 - 36 6i 2 2 - 171 a 2 i& 02 2 
+ 5 1 3 6 30 6 2 2 + 537 a u % 2 2 + 351 b 02 2 b 12 - 271 a u b 02 3 + 542 b 20 b 02 3 
+ 756 b 20 b 30 b 02 + 2268 &20&30O11 - 20 b 02 4 + 1120 6 20 4 + 798 a u 2 b 20 2 

- 2240an6 2 o 3 - 1512 6 2 o0 4 o + 1008 6 2 o 2 o 2 i - 252 b 20 2 b 12 
+ 1806 b 20 % 2 2 + 2240 b 20 3 b 02 



Sec.10.2). 



x = -y + a n xy + . . . + a n _ 1A x n 1 y 

y =x + b 20 x 2 + b 02 y 2 + ... + b n . 2>2 x n ~ 2 y 2 + b nfi x 




(33) 
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To solve the systems of polynomials equations we use Grobner bases. Solving system of 
9 (non-zero) equations of Sys(9) requires higher performance computers and ours are not 
up to it. 

A careful analysis of polynomials from Sys(m) shows that for any m they always 
are weighted-homogeneous. For example, the polynomial P 2 given by (pM|) is weighted- 
homogeneous if we give weight 2 for a 2 i, &12 and 6 30 , and weight 1 for the remaining 
variables. 

We observe that all polynomials from Sys(m) are weighted- homogeneous if we choose 
the following weights 

1 . i + j — 1 for parameters and fry 

2. 2i + 1 for c 2 i+i- 

We introduce new parameters B^, and C 2 j+i putting 

A % ij 3 1 = a>ij B l ^ j 1 = b^, C%Hl = c 2i+1 (35) 

After this reparametrization system (133]) reads 

x = -y + A n xy + ... + A^x^y 
y = x + B 20 x 2 + B 02 y 2 + ... + B^x^y 2 + B^x 

As in the case of isochronous center the Urabe function is odd, we search it under the form 

00 

h(X) = ClH\X 2k+l = dX + C|A 3 + C 5 5 A 5 + C 7 7 X 7 + ... 

k=0 

By a simply use of (|35|) . from the isochronicity conditions for system (136|) . expressed in 
terms of its parameters {A?} and {B rs }, it is easy to recover the parameters values {a^} 
and {b rs } when system (1331) admits isochronous centers at the origin O. 

The described reparametrization gives rise to homogeneous equations and reduces the 
number of parameters appearing in fl36l) by one. First we assume B 2 q = 0, and we solve 
the isochronicity problem for system fl36l) under this assumption. Next, for B 2 q 7^ 0, we 
apply to system fl36|) the following change of coordinates 

(z,V) | - > (-5-1 -5-) 
-D20 -D20 
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Where 

> 

Hence, without loss of generality we can put B 2 o = 1, and find the parameters values for 
which the center is isochronous. 

Note that the resolution of the polynomial system issued from the 19 derivations and 
associated eliminations for system (|15|) (with 9 parameters), exceed our computer facilities. 

7.3 On efficiency of C-algorithm 

Usually for the search of isochronous centers one uses the method of normal form (see for 
exemple [XT] [7]). Thus it is interesting to compare C-algorithm with it, when there are 
applied to systems reducible to Lienard type equations. 

As an example for such investigations we choose the Abel systems (1251) with 2 < n < 9. 
The normal form (NF) algorithm used is the one described in [20] which is universal and 
efficient. As explained in point (3) of Sec. 7.1 for the system fT25|) with n parameters one 
computes coefficients Ci, c 3 , . . . , c 2m _i of the Urabe function when C-algorithm is used. 
When the normal form method is used one computes its first 2n + 1 terms. 

The results are presented in the table where the time unit is one second. 
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Table 1: CPU time on Pentium 1,46 GHz 



n 


C- algorithm 


NF algorithm 


2 


~ 


0,060 


3 


0,001 


0,160 


4 


0,004 


0,784 


5 


0,008 


4,728 


6 


0,016 


31,430 


7 


0,052 


263,033 


8 


0,116 


2335,962 


9 


0,284 





The superiority of C-algorithm is obvious. 
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